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Abstract 

Any pure three-qubit state is uniquely characterized by one phase and four positive parameters. 
The geometric measure of entanglement as a function of state parameters can have different ex- 
pressions. Each of expressions has its own applicable domain and thus the whole state parameter 
space is divided into subspaces that are ranges of definition for corresponding expressions. The 
purpose of this paper is to examine the applicable domains for the most general qubit-interchange 
symmetric three-qubit states. First, we compute the eigenvalues of the non- linear eigenvalue equa- 
tions and the nearest separable states for the permutation invariant three-qubit states with a fixed 
phase. Next, we compute the geometric entanglement measure, deduce the boundaries of all sub- 
spaces, and find allocations of highly and slightly entangled states. It is shown that there are three 
applicable domains when the phase factor is n/ 2 while other cases have only two domains. The 
emergence of the three domains is due to the appearance of the additional W-state. We show that 
most of highly entangled states reside near the boundaries of the domains and states located far 
from the boundaries become less-entangled and eventually go to the product states. The neighbors 
of W-state are generally more entangled than the neighbors of Greenberger-Horne-Zeilinger(GHZ) 
state from the aspect of the geometric measure. However, the range of the GHZ-neighbors is much 
more wider than the range of the W-neighbors. 
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I. INTRODUCTION 



Entanglement is a property of quantum states that does not exist classically. Two or 
more subsystems of a quantum system are said to be entangled if the state of the entire 
system cannot be described in terms of a state for each of the subsystems \V]. This property 



of composite quantum systems, which exhibits quantum corre 
a resource for many processes in quantum information theory 



ations between subsystems, is 



3|, |j, |5| . Since theprofound 



measures of entanglement, i.e. the entanglement of formation and distillation 



flaaa 



have 



not been properly generalized to multipartic 
entanglement via other measures [l(J 



11 



12 



e systems, the study of quantifying multipartite 
13L Q is a necessity. 



The entanglement of a given pure state can be characterized by a distance to the near- 
est unentangled state 15|. A whole class of such entanglement monotones, based on the 



Euclidean distance of a given multipartite state to the nearest fully separable state, was 
constructed in Ref.flq]. Subsequently, a geometrically motivated measure of entanglement, 
known as geometric measure, was introduced by Wei and Goldbart [jjj. It is a decreasing 
function of the maximal overlap P max and is suitable for any partite system regardless of its 
dimensions. The maximal overlap has several different names and we list all of them for the 
completeness: maximal probability of success [13| . entanglement eigenvalue [l7|, infective 
tensor norm the largest Schmidt coefficient 19] and maximum singular value 20]. 
The geometric measure has an advantage that it can be computed analytically for multi- 



parameter states. Recently, explicit expressions for the maxima! 



for three- 



2d, 



21 



22 



231 ] as well as for multi-qubit states 



24 



overlap have been derived 



25 



26 



271 ]. It turned out 



that the maximal overlap, depending on coefficients of a quantum state in a computational 
basis, can take two different values. It is equal to either the square of the largest coefficient 
or the square of the circumradius of a cyclic polygon constructed by the coefficients of the 
quantum state. This means that the whole parameter space is divided into two subspaces 
each of which has its own expression for the geometric measure. 

In spite of these achievements, still we lack sufficient knowledge to classify generic three- 
qubit pure states by the geometric measure. They have five local unitary(LU) invariants 



including four positive parameters and a gauge phase 7 [19|, [28|, |29| . The maximal overlap 
of these states is not known yet. Only three-qubit states which are expressed as linear 
combinations of four(or less) orthogonal product states have been considered so far 22]. 
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In fact, all of these states have real coefficients because the phases of their coefficients can 
be eliminated by LU-transformations. Thus, the contribution of the gauge phase to the 
maximal overlap has remained a mystery. On the other hand, the most recent results 30] 
have shown that the gauge phase plays an important role. It parameterizes the family of 
maximally entangled states and identifies W-class pure states with the boundary of pure 
states. 

In this paper we would like to take into complete account the effect of the gauge phase 
in the geometric measure of entanglement. We compute the maximal overlap as well as the 
nearest product states for a given value of the gauge phase. We will show in the following 
that depending on the phase factor 7 the whole parameter space is divided into the two 
or three domains, each of which has a particular expression for the geometric measure. In 
addition, we will show that most of highly entangled states reside near the boundaries of the 
domains. We will call these highly entangled states as GHZ-neighbors. The states located 
far from the boundaries become less-entangled and eventually go to the product states. But 
there is different kind of the highly entangled states. These states reside around W-states. 
We will call these highly entangled states as W-neighbors. The W-neighbors are generally 
more entangled than the GHZ-neighbors from the aspect of the geometric measure. However, 
the range of the GHZ neighbors is much more wider than the range of the W-neighbors. 

The paper is organized as follows. In section II following Ref. j2l| we transform the 
nonlinear eigenvalue equations into the Lagrange multiplier equations. In section III we 
solve the Lagrange multiplier equations analytically for 7 = and 7 = tt/2. It turns out 
that both cases give five different eigenvalues. Also every eigenvalue has its own available 
region in the parameter space. In section IV we compute the geometric measure for 7 = 
case. It turns out that two of the five eigenvalues contribute to the geometric measure. 
This means that the whole parameter space is divided into two applicable domains. In 
section V we compute the geometric measure for 7 = tt/2 case. It is shown that the whole 
parameter space is divided into the three applicable domains. In section VI we compute 
the eigenvalues and the geometric measure for 7 = 7r/4 numerically. It is shown that when 
7 = 7r/4, there are six different eigenvalues. However, only two eigenvalues contribute to the 
geometric measure. In section VI a brief conclusion is given. In appendix we have shown that 
Lagrange multiplier equations for arbitrary 7 provides a solution whose multiplier constant 
is zero. 



3 



II. GENERAL FORMALISM 



In this section we clarify our notations, give necessary definitions, define three-qubit sym- 
metric states and transform nonlinear stationarity equations to a system of linear equations. 

A. Preliminaries 

The maximal overlap of n-qubit pure states is given by 

Pmax = max |((gi|0(g2|0- • •0(gn|)|V , )| 2 , ( 2 -l) 
91,92,— ,9n 

where the maximization is performed over single qubit pure states. Constituents \qi), |<?2), 
\q n ), the nearest product state from can be computed via the non-linear eigenvalue 
equations 

• ■■(in-i\^) = tk\q n ), (qi\ ■ ■ ■ (qn-2\(q n \tp) = tii\q n -i), (Q2\ ■ ■ ■ (q n \*P} =m\qi), 

(2.2) 

where /ij's are the eigenvalues of Eq. (12.21) . Then the geometric measure G of the quantum 
state is defined as G(ip) = 1 — P maxi where P ma x = max(/i-). 
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For simplicity, we take a quantum states which possess a permutational symmetry 



31 



331 ]. These states have three independent parameters and, through an appropriate LU 



transformations, can be brought into the symmetric form [1{ 

= g\Q00) + t|oil) + t|101) + t|110) + e i7 /i|lll), (2.3) 



where we follow the notation of Ref. [30(. In above equation all coefficients g, h and t are 
positive and satisfy the normalization condition g 2 + 3t 2 + h 2 = 1. The phase 7 has the period 
7r and ranges within the interval —tt/2 < 7 < n/2. Note that Eq. (12.31) is not a Schmidt 
decomposition for l^) since the Schmidt normal form imposes additional conditions(namely, 
a lower bound on g) on state parameters. We would like to abandon these additional 



constraints and apply the general method proposed in Ref. 



21] to symmetric states Eq. (12.31) . 



B. Modified stationarity equations 

In this subsection we would like to present the method for solving stationarity equations 
for the quantum state given in Eq. (12.31) . In the case of three-qubit pure states the method 



developed in Ref. 2l| transforms the system of nonlinear equations to a system of linear 
equations. In spite of this essential simplification, it is impossible to get analytic expressions 
for generic three-qubit states since the solution of the linear eigenvalue equations reduces 



to the root finding for a couple of algebraic equations of degree six 22]. However, the 
permutation symmetry of reduces this pair of algebraic equations to a single algebraic 
equation of degree six. Furthermore, there is a solution which holds for all values of state 
parameters 30|. The separation of this global solution allows us to solve explicitly the 
eigenvalue equations for 7 = and 7 = it/2 and leads us to a quartic equation for remaining 
cases. The quartic is the highest order polynomial equation that can be solved by radicals in 
the general case. But expressions for roots are impractical and we will carry out numerical 
analysis instead. 

The method enables us to express eigenvalues p 2 via the reduced densities p A , p B and 
pAB Q f qubitg A and B in a form: 



p = - max (1 + r 1 ■ si + r 2 ■ s 2 + Gi j s li s 2 j] 

4 |si| = |s 2 |=l 



where 



n = Tr (pV) , r 2 = Tr (p B a) , G l3 = Tr (p^a, ® a 3 ) 
and crj's are Pauli matrices. Explicit calculation shows 



(2.4) 



(2.5) 



r = n = r 2 = (2ht cos 7, 2ht sin 7, g-h- tr) (2.6) 

f 2t(g + t) -2/itcos7 X 

Gij = -2t(g-t) -2/itsin7 

^— 2/itcos7 — 2/itsin7 g 2 + h 2 — t 2 J 

It is worthwhile noting that r\ is identical with r 2 and is a symmetric matrix. These 
properties arise due to the fact that we have chosen the symmetric state in Eq. fl2.3p under 
the qubit-exchange. As will be shown in the following these properties drastically simplify 
the calculation procedure. Since r%, r 2 and G^ are explicitly derived, the eigenvalues p 2 can 



be computed if Si and s 2 are known. Due to the maximization in Eq.f l2.4p these vectors can 
be computed by solving the Lagrange multiplier equations: 



ri + Gs 2 = X 1 s 1 



r 2 + G T s 1 = \ 2 s 2 



(2.7) 
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where the superscript T stands for transpose and A;'s are the Lagrange multiplier constants. 
From the properties r% = r 2 and Gij = Gji Eq. (12.71) can be reduced to a single equation 



r + Gs = \s (2.8) 
where A = Ai = A2 and s = si = S2. Letting 

s = (sin 9 cos ip, sin 9 sin ip, cos 9), (2-9) 

Eq. (12.81) reduces to 

2ht cos 7 + 2t(g + t) sin 9 cos <p — 2ht cos 7 cos 9 = A sin 9 cos <p (2.10a) 

2ht sin 7 — 2t(g — t) sin 9 sin ip — 2/it sin 7 cos 9 = X sin # sin 9? (2. 10b) 



(g 2 — t 2 )(l + cos#) — h 2 (l — cos9) — 2/itcos7sin^cosy9 — 2 /it sin 7 sin 6 1 sin = Acos#. (2.10c) 

Solving 9, ip and A from Eq. (l2.10|) . one can compute the eigenvalues for the symmetric 
canonical state (12.31) by inserting the solutions into Eq. (12.41) . In the next section we will 
solve analytically Eq. (l2.10p at the particular phases 7 = and 7 = n/2. By making use 
of the solutions we will compute fii and P ma x = max(/xf) for the corresponding quantum 
states. 

III. EIGENVALUES 

In this section Eq. (12. 101) will be solved at 7 = and tt/2 separately. Since numerical 
calculation is needed to analyze the 7 = 7r/4 case, we deal with this case in different section 
(see section VI). 

A. 7 = case 



For this case Eq. (12.101) reduces to 

2t(g + t) sin^cos(/9 + 2ht(l — cos6*) = A sin 9 cos <p (3.1a) 

— 2t(g — t) sin 9 sin p = A sin 9 sin <p (3.1b) 

(g 2 -t 2 )(l + cos9) - h 2 {l - cos9) - 2htsin9cosp = \cos9. (3.1c) 
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Eg. (13. lbf ) implies that the solutions for the 7 = case are categorized by 9 = 0, if = 0, 
if = it and A = —2t(g —t) 1 . 



1. 8 = case 



When 9 = 0, Eq. fl3.lap and Eq. (13. lbj) are automatically solved, and Eq. fj3.lcl) gives 



A = 2(^-0- (3-2) 

Now s = (0, 0, 1) and Eq. (12.41) together with the normalization condition g 2 + 3t 2 + h 2 = 1 
gives the eigenvalue 

/4 = g 2 - (3.3) 

2. if = case 

For this case Eq. ( 13. lbj) is automatically solved and the remaining equations are 

2t(g + t) sin 6 + 2ht{\ - cos 9) = A sin (3.4a) 

(c? 2 - t 2 )(l + cos^) - h 2 (l - cos 9) - 2htsin9 = \cos9. (3.4b) 
Since sin(0/2) ^ 0, Eq.fJH3a]) reduces to 

A = 2htz + 2t 2 + 2tg (3.5) 

where z = tan((9/2). Inserting Eq. fl3.5p into Eq. fl3.4bj) . one can derive an equation 



(hz + g + t) (tz 2 - hz + g - 2t) = 0. (3.6) 



Eq. (13.61) implies that the ip = case is also categorized again by following three cases: 

g + 1 r + r_ 



(3.7) 



h ' 2t' 2t 
where 

r± = /i ± v / /i 2 + 4t(2t-^). (3.8) 



1 The case 8 = tt can be excluded by Eq. (l3.1al) . 
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First, let us consider the case of z = 



— (g + t)/h. In this case Eg. (13.51) gives 



A = 0. (3.9) 

Since, in this case, 



. 2h(g + t) h 2 -{g + tf 

s * = sm6 = - h> + ( 9+ tr Sy = 0, Sz = h*+( g +tr (3 0) 



it is straightforward to compute the eigenvalues for this case, which is 

_ g 2 h 2 + t 2 (g + t) 2 

h 2 + ( g + t )2 ■ V- LL ) 

Next, let us consider the case of z — r±/2t simultaneously. In these cases Eq. (13.50 gives 

A = hr± + 2t(g + t). (3.12) 

Since, in these cases, 

4tr ± r| - At 2 

s y = 0, s z = - (3.13) 



x ri + At 2 ' y ' z r 2 ± + At 2 

one can show directly that the eigenvalues are 

_ (hr ± + At 2 ) 2 

Since z = tan(^/2) should be real, the eigenvalues ia\ are available only when 

h 2 

g<2t+—. (3.15) 

At v ' 

3. ip = 7r case 

For this case Eq. f!3.1bl) is automatically solved and the remaining equations are 

- 2t(g + t) sin 6 + 2ht(l - cos B) = -A sin 9 (3.16a) 

(g 2 - t 2 )(l + cos 9) - h 2 (l - cos 9) + 2htsin9 = \cos9. (3.16b) 

Since Eq. (l3.16|) can be derived from Eq. (13.41) by changing 9 — ► —9, the solutions for this 
case are also categorized by 

g + 1 r + r_ 
Z ~ h ' ~2t' ~2t 
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(3.17) 



Since Eq. (l3.16al) reduces to 

A = -2htz + 2t 2 + 2tg, (3.18) 

comparison of Eq. (l3.18p with Eg. (13. 5 p shows that the Lagrange multiplier constant A is same 
with the case of <p = 0. Since, furthermore, s x = sin 9 cosy? and s z = cos 9 are invariant 
under 9 — > —9 and <p = — > 9? = 7r, this fact implies that the eigenvalues for this case are 
exactly same with those for p = case. 



4- A = 2t 2 — 2gt case 

For this case Eq. (13.1bl) is automatically solved and the remaining equations are 

2t(g + t) sin 9 cos <p + 2ht{\ - cos 9) = -2t(g - t) sin 9 cos cp (3.19a) 

(g 2 -t 2 )(l + cos9) -h 2 (l -cos9) - 2ht sin 9 cos <p = -2t(g - t) cos0. (3.19b) 
Since Eq. fl3.19al) gives a relation 

hi — cos 6* 

cos(^ = -- — — , 3.20 

2g sin 6* 



combining Eq. (I3.19bl) and Eq. (13.201) enables us to express cos 9 and sin 9 as 



£ 2 -fr 2 + fft . , ^Ag(g + 2t){V + gt) 
cos9 = — smS = ±- — . 3.21 

g 2 + h 2 + 3gt g 2 + h 2 + 3gt K ' 

For a time being we choose the upper sign in sin^. Then, Eq. (l3.20p reduces to 



h / q + 2t , 

cosp = -J— -. 3.22 

2\jg(h 2 + gt) K } 

At this stage it is worthwhile noting that the eigenvalue in this case is available when 

(3g - 2t)h 2 + AgH > (3.23) 
because of —1 < cosp < 1. Of course, the corresponding sinp is 



smv=± V • (3 - 24) 

Again we choose the upper sign in smp. Then, it is straightforward to compute s, whose 
components are 

h(g + 2t) y/(g + 2t)(3gh 2 + Ag 2 t - 2h 2 t) g 2 - h 2 + gt 

Sx ~ ~g 2 + h 2 + 3gt Sy ~ g 2 + h 2 + 3gt §z ~ ~ g 2 + h 2 + 3gt 

(3.25) 
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Inserting Eq. fj3.25j) into Eq. (12.41) gives the eigenvalue for this case as follows: 

g{gh 2 + At 3 ) 



(3.26) 



g 2 + h 2 + 3gt' 

It is easy to show that the choice of other sign in sin# and sinp does not change the 
eigenvalue /i^. 

The eigenvalues for 7 = case are summarized in Table I. 



name 


eigenvalue 


A 


available region 




9 2 


2{g 2 -t 2 ) 


all 


n\ 


g 2 h 2 +t 2 (g+t) 2 
h 2 +(g+t) 2 





all 


& 


(hr + +At 2 ) 2 
i :2 _+At 2 


hr+ + 2t(g + t) 


g<2t + h 2 /(At) 




(hr-+At 2 ) 2 
r 2_ +4t 2 


hr„ + 2t(g + t) 


g<2t + h 2 /(4t) 


i4 


g{gh 2 +U 3 ) 
g 2 +h 2 +3gt 


2t{t - g) 


(3g - 2t)h 2 + AgH > 



Table I: Eigenvalues for 7 = case 



B. 7 = 7r/2 case 

For this case Eq. (12.101) reduces to 

2t(g + t) sin 6 cos <p — A sin 9 cos <p, (3.27a) 

— 2t(g — t) sin 9 sin ip + 2ht(l — cos 9) = A sin 9 sin ip, (3.27b) 

(g 2 -t 2 )(l + cos9) - h 2 (l — cos 9) — 2ht sin 9 sin <p = A cos 9. (3.27c) 

Eq. fl3.27ap guarantees that the solutions for this case are categorized by 9 = 0, ip = 7r/2, 
<p = 3tt/2 and A = 2t(g + 1). Since the calculation procedure for the first three cases are 
similar to the 7 = case, we will briefly sketch the final result only. Although the calculation 
procedure for the last case is also similar to the previous case, it gives a non-trivial available 
region, which is important to compute the geometric measures in next section. Therefore, 
we will present the last case in detail. 

When 9 = 0, the Lagrangian multiplier constant is same with Eq. fl3.2p and the corre- 
sponding eigenvalue is 

v 2 P = g 2 . (3.28) 
10 



When p> = 7r/2, there are three types of solutions depending on z = tan(0/2). If z = (g—t) / h, 
we have vanishing Lagrange multiplier constant and the corresponding eigenvalue is 

_ g 2 h 2 + t 2 (g-t) 2 

Vl ~ h * + (g-ty • (3 ' 29) 

When z = s±/2t, where 

s± = h± ^h 2 + 4t(2t + g), (3.30) 

the corresponding Lagrange multiplier constants are hs± — 2t(g — £), and the corresponding 
eigenvalues are 

_ (hs ± + At 2 ) 2 , 31) 

It should be noted that v\ are available in entire parameter space, while fi± in 7 = case 
is restricted by Eq. fl3.15l) . As in the case of 7 = 0, ip = 3n/2 case does not give a new 
eigenvalue. This case just reproduces v\ and v±. 

Finally, let us discuss A = 2t(g + t) case. For this case Eq. fl3.27al) is automatically solved 
and the remaining equations are 

2ht(l — cos 6) — 2t(g — t) sin 9 sin <p = 2t(g + t) sin 9 sin <p (3.32a) 

{g 2 - t 2 )(l + cos6) - h 2 (l - cos9) - 2ht sin 6 sin p = 2t(g + t) cos 6. (3.32b) 
Since Eq. (l3.32al) gives a relation 

hi — cos 9 

sm<p = - — — , 3.33 

2g sin 9 

combining Eq. (l3.32bl) and Eq. (13.331) yields 

g 2 -h 2 - gt 

COSe = - gi + h*-3gt (3 ' 34) 
The requirement — 1 < cos 6 < 1 gives first available condition 

(g-2t)(h 2 -gt) > 0. (3.35) 

Now we choose sin 9 as 

. ^Ag{g-2t){W-gt) 

Sm6= g 2 + hi-3gt • (3 - 36) 
Then from Eq. fl3.33p sinp becomes 



h / g-2t 
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Another requirement — 1 < simp < 1 gives second available condition 

(g - 2t){3gh 2 - Ag 2 t + 2h 2 t) > 0. 

Choosing cosy? as 



cosy? 



'3gh 2 - AgH + 2hH 



Ag(h 2 -gt) 

it is straightforward to show that the eigenvalues for this case is 

2 g(gh 2 -4t 3 ) 



(3.38) 



(3.39) 



(3.40) 



2 g 2 + h 2 -3gf 

It is easy to show that the different choices in the sign of sin 8 and/or cosy? do not change 
the eigenvalue. Although the available region for v\ is restricted by Eq. (13.351) and Eq. (13.381) . 
one can show that Eq. ( 13.381) implies Eq. ( 13.351) already. To show this explicitly let us consider 
g > 2t case first. In this case Eq. (13.381) imposes h 2 > 4g 2 t/(3g + 2t). Therefore 

V* , gt 



h - gt> 



gt 



-(g-2t)>0. 



3g + 2t ° 3g + 2t 

Similarly, one can show that Eq. (13.381) implies Eq. fl3.35p for g < 2t region too. Therefore, 
the available region for v\ is restricted by Eq. ( 13.381) only. 
The eigenvalues in 7 = ir/2 case is summarized in Table II. 



name 


eigenvalue 


A 


available region 


vl 


g 2 


2{g 2 - t 2 ) 


all 




g 2 h 2 +t\g-t) 2 





all 


-I 


(hs + +At 2 ) 2 
s\+At 2 


hs + - 2t(g - t) 


aU 


v 2 _ 


(hs-+4t 2 ) 2 
s 2 _+At 2 


hs^ - 2t(g - t) 


all 


v\ 


g(gh 2 -4t 3 ) 
g 2 +h 2 -3gt 


2t{g + 1) 


{g - 2t)(3gh 2 - AgH + 2hH) > 



Table II: Eigenvalues for 7 = ir/2 case 



C. h -> limit 

Since is independent of 7 in the /z — ^ limit, all eigenvalues for 7 = and 7 = 71/ 2 
cases should be same including the available region in the parameter space. Note that 
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= and v\ = v 2 _ in the h — > limit. In this limit the eigenvalues for 7 = exactly 
coincide with eigenvalues for 7 = 7r/2 as following: 

4+3 4 |3 

2 2 222. 2 22 ™ 2 2 ^ L /o At\ 

fIp = Up = g fj ll = Ul =t n 2 = v ± = — — n ± = v 2 = —— -. (3.41) 

In addition, first three eigenvalues in Eq. (13.411) are available in the full parameter space 
and the last one is available only at g < It. Thus, our calculational results are perfectly 
consistent in the h — > limit. 



IV. GEOMETRIC MEASURE FOR 7 = 

In this section we would like to compute the geometric entanglement measure defined 

G(V0 = l-P max W (4-1) 

for 7 = case. In order to compute P ma x we would like to emphasize three points, which 
simplify the following calculation. Firstly, note that P ma x is given by 

P max = max(X 2 ). (4.2) 

Therefore, we should choose the largest eigenvalue from all eigenvalues, each of which has 
its own available regions in the parameter space. Secondly, note that 

2 2 128/it 7 / 2 / h 2 \ 3/2 

wrmwr^) { 2t+ T t - g ) • (4 - 3) 

This means that /x 2 , is always smaller than f£ in the available region g < It + h 2 /(4t). 
Therefore, we can exclude ix 2 _ from beginning for the computation of P m ax- Thirdly, note that 
Pmax is obtained from the eigenvalues whose Lagrange multiplier constants are positive 2l| . 
This fact excludes /j, 2 too. Considering all of these facts and available regions, it is convenient 
to divide the whole parameter space into the following four regions: 

h 2 

(region I) g > 2t + — : P max = /4 ( 4 - 4 ) 

h 2 

(region II) t < g < 2t + — : P max = max(/i P , 
(region III) g < t & C x > : P max = max(/i^ , 
(region IV) g < t h C x < : P max = /2 2 + 
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V 



FIG. 1: (Color online) Fig. la is a plot of the applicable domains in (u, u)-plane for 7 = 0. The 



this domain is around large g region. Fig. lb is plot of (u, v )-dependence of P ma x for 7 = 
case. Many highly entangled states are represented as a valley in this figure. Around u ~ and 
(u ~ ir/2,v ~ 0) there are a lot of less entangled states. To compare the applicable domains with 
Pmax w e plot both simultaneously in the (u, v) plane in Fig. lc. The black thick line is a boundary 
between domains. The blue-color and white-color represent the highly- and less-entangled states 
respectively. Fig. lc shows that the highly-entangled states reside around the boundary between 
domains. 



principal domain P, 



max — 



fj? p is located in small v and large u region. This fact indicates that 



where 



d = (30 - 2t)h 2 + AgH. 



(4.5) 
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In order to compare ji 2 + with /if, we compute fi+ — which is 

2 



2 2 
fl+-fi 2 



{r\ + 4t 2 )(g 2 + h 2 + 3gt) 



(a! + (3^h 2 + 4t(2t- g) 



where 



(4.6) 



(4.7) 



ai = h 6 + ghH + 8hH 2 + 20gh 2 t 3 + lQgH 4 + Ah 2 t 2 (2t 2 - g 2 ) 
fa = h(h 4 + 3gh 2 t + Agh 2 + Ahh 2 + 8gt 3 ). 

Since the last term in a%, 4h 2 t 2 (2t 2 — g 2 ), is non-negative in the region g < t, both a\ and 
(3\ are non-negative in region III. In region III, therefore, P m ax becomes 

In region II it has been shown in Ref. 30( that /i|> = when T>i = 0, where 

V 1 = gh 2 -(g + t) 2 (g-2t). (4.8) 

Therefore, the region II should be divided into two regions, i.e. T>\ > and T>\ < 0. Simple 
consideration shows that \j? p > n\ when T>\ < and //|> < fi 2 ^ when Pi > 0. Combining all 
of these facts, one can conclude 



h 2 

(region A) 5 > 2t + — : P n 



Hp 



(4.9) 



(region B) t < g < 2t + — & P a < : 
(region C) t < g < 2t + — & V 1 > : 



p 



p 

J. 77 



Hp 



M 2 + 



(region D) g <t : P„ 



Now, we would like to unify the regions as many as possible to simplify the expression 
of Pmax- First, one can show that T>\ is always non-positive in region A as following. Since 
h 2 < 4t(g — 2t) in region A, in this region 

Pi = gh 2 -(g + t) 2 (g - 2t) < -(g - 2t)(g - t) 2 < 0. (4.10) 

Second, one can show easily that D\ is always non-negative at region D as following. In this 
region 

V l =gh 2 + (g + t) 2 (2t-g)>0 (4.11) 

because both terms are non- negative. Combining these facts and Eq. fl4.9p makes P ma x to be 
expressed as 

Hp when T>i < 
h\ when T> x > 0. 



(4.12) 
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In order to understand the behavior of P ma x more clearly we introduce the two parameters 
u and v as following: 



with < u,v < tt/2. Then, one can plot the applicable domains T>\ < and T>\ > in the 
u — v plane, which is Fig. la. As Fig. la has shown, the domain for T>\ < is biased in 
the small v and large u region. This indicates that the domains for T>\ < is around large 
g region. The remaining region is the domain for T>\ > 0. As will be shown in next section, 
the number of the applicable domains for 7 = 71/ 2 case is not two but three. This means 
that the phase factor 7 has great impact in the geometric measure of entanglement. 

Fig. lb is (u, v )-dependence of P ma x given in Eq. fl4.12l) . At u — 0, which means h = 1, 
Pmax becomes 1 because it is separable state. At v = and u = 7r/2, which means that 
g = 1, Pmax becomes 1 again. Between them there is valley, which represents the set of 
the highly entangled states. There is different kind of the highly entangled states around 
u = v = tt/2. These highly entangled states are states located near W-state, \W) = 



In order to compare P max with the applicable domains we plot P ma x and the boundary 
of domains simultaneously in u — v plane in Fig. lc. In Fig. lc the black thick line is a 
boundary of the domains. The thick-blue color and light-blue (or white) colors represent 
the highly-entangled and less-entangled states, respectively. In the right-upper corner there 
are many highly entangled states which are located near W-state. Another type of the 
highly entangled states reside near the boundary of the applicable domains. Apart from the 
boundary more and more the quantum states lose the entanglement, and eventually reduce 
to the separable state. 

Now, we consider several special cases. First example is t = l/y/3 and g = h = 0. In 
this case T>i = 2^3/9 > and r+ = ^/8/3, which gives Pmax = 4/9. Second example is 
t = and g > h. In this case T>\ = —g(g 2 — h 2 ) < and P m ax = g 1 ■ Third example is 
t = and g < h. In this case T>\ = g(h 2 — g 2 ) > and r + = 2h, which gives P max — h 2 . 
The second and third examples are consistent with P max (GH Z) = max(|a| 2 , |/3| 2 ), where 
\GHZ) = a|000) + /3|111). Fourth example is g = case. In this case V x = 2t 3 > and 
r + = h + \/h 2 + 8t 2 , which results in 



g = sinucosw, 



t = sin u sin vj h = cos u 



(4.13) 



(1/V3)(|011) + |101) + |110». 



P 



(h 4 + 8h 2 t 2 + 8t 4 ) + h(h 2 + 4t 2 )Vh 2 + 8t 2 



(4.14) 



max 



(h 2 + Qt 2 ) + hVh 2 + 8t 2 
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One can show that various limits of Eq. fl4.14p are consistent with the previously derived 
results. The last example is h = case. In this case it is easy to show 

{q 2 when q > 2t 

(4.15) 
4t 3 /(3t - g) when g < 2t. 



Eg. (14.151) is perfectly in agreement with the result of Ref. 22]. 



V. GEOMETRIC MEASURE FOR 7 = vr/2 

In this section we would like to compute the geometric entanglement measure for 7 = 71/ 2 
case. From the constraint of the positive Lagrange multiplier constant we can exclude v\ 
and v 2 _ from beginning stage for the computation of the geometric measure. Next, we should 
examine the sign of the Lagrange multiplier constant for z/i, that is 

A + = hs+-2t(g-t). (5.1) 

It is easy to show that A + > in g < t region. Also it is straightforward to show that 
A + > when C + > and A + < when C + < 0, where 

C + = h 2 (2g + t)-t(g-t) 2 . (5.2) 
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c + < o v C+ > 



c 2 < o 



- y C 2 > 

K ^ 



FIG. 2: Pictorial representation for C2 > 0, C 2 < 0, C+ > 0, and C+ < when g >t. 

Examining Table II and Eq.f l5.2p leads us to divide the whole parameter space into the 
following ten regions: 

(z) g>2t (5.3) 



where 





$ion 


I) C 2 < & C+ < : 


P — 
1 max 




ref 


$ion 


II) 


C 2 > & C+ < : 


P 

1 max 


- max (i/p, U2) 


rej 


lion 


HI) 


C 2 < & C+ > : 


P 

1 max 


= max(i/p, v+) 


rej 


lion 


IV) 


C 2 > & C+ > : 


P 

1 max 


= max(z/p, u^., v 








(m) t<g<2t 






rej 


j;ion 


V) 


C 2 > & C+ < : 


P 

1 max 


-4 


ref 


$ion 


VI) 


C 2 < & C+ < : 


P 

1 max 


= max(z/p, v%) 


rej 


$ion 


VII) 


C 2 > & C+ > : 


P 

1 max 


= max(i/p, u.) 


rej 


$ion 


VIII) 


C 2 < & C+ > 


■ Pmax — max(z/p, v + , 








(iii) <? < t 






ref 


$ion 


IX) 


C 2 <0: P max = max(z4, v\) 




ref 


$ion 


X) 


C 2 > : P max — v\ 







C 2 = (3g + 2t)h 2 - AgH. (5.4) 



Although the whole space is divided into the ten regions, one can show that some regions 
do not exist. In order to show this it is convenient to introduce 

3g + 2t + 2g + t K ' 
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c+ < 



c 3 < o 



c 2 < o 



h + /i 3 h 2 
t<g<2t 



C + > 
C 3 > 

c 2 > o 

h 2 



C 4- < 



c 2 < o 



c 3 <o 



^+ ^2 ^3 



C+ > 

c 2 > o 
c 3 >o 

+ /J 2 



FIG. 3: Pictorial representation for C 2 > 0, C 2 < 0, C+ > 0, C+ < 0, C 3 > and C 3 < when 
t<g<2t (Fig. 2 a) and g > 2t (Fig. 2 b). 



Then, their difference becomes 



ho - h. 



t(g + tf 



-(5g-2t). 



(5.6) 



(3g + 2t)(2g + t) 

Eq. (15.61) implies that h 2 > h + in the region g >t. Then the regions C 2 > 0, C 2 < 0, C + > 0, 
and C + < when g > t can be represented as Fig. 2. With an help of Fig. 2 it is easy to 
understand that there is no region which satisfies both C 2 > and C + < when g >t. This 
implies that region II and region V do not exist in the whole parameter space. 
In order to compare up with v\ we compute vp — u 2 , which is 



g{g + t){g-2tf 



(5.7) 



* 1 g 2 + h 2 -3gt ' 
Therefore, the sign of v\ — v\ is determined by g 2 + h 2 — 3gt. If C 2 > 0, h 2 > h 2 and 

9 Vtf-3 9t >M9_^|±l (5.8) 

Therefore, if C 2 > in g > 2t region, vp > z/f. Thus, we can exclude z/f in region IV. 
Similarly, one can show that if C 2 < in t < g < 2t region, v 2 p < z/f. Therefore, we can 
exclude vp in regions VI and VIII. 



Next, we compute v P — v+, which is 

2 



Z/p - V, 



s 2 + At 2 



(a 2 + (3 2 ^/h 2 + 4t{2t + gfj (5.9) 
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where 



a 2 = -h 4 + (g + 2t)(g - U)h 2 + 2t(g - t)(g + 2t) 2 (5.10) 
(3 2 = h(g 2 -h 2 -At 2 ). 

Direct calculation shows that in g > t region up = u\ when C 3 = 0, where 

C 3 = gh 2 -(g-t) 2 (g + 2t). (5.11) 

In addition, simple consideration shows that in g > t region up > u^_ when C3 < and 
v 2 p < u 2 when C 3 > 0. 

In order to check which eigenvalue is dominant in each region it is convenient to introduce 
another parameter 

k3=(g -tn 9+ 2 t ) 

9 

Then, it is easy to show 

h+<h 2 <h 3 when 2t < g (5.13) 

h + < h 3 < h 2 when t < g < 2t. 

Eq. fl5TT3l enables us to represent C 2 > 0, C 2 < 0, C+ > 0, C + < 0, C 3 > and C 3 < 
in one-dimensional coordinate, which is illustrated in Fig. 3. With an help of Fig. 3 one 
can show easily that in region III C 3 is always non-positive and therefore, P ma x becomes up. 
Using Fig. 3a P max in region VII is u\. Using Fig. 3b again one can show that region IV is 
divided into 

(region IV-a) C 2 > & C 3 < : P max = v 2 P (5.14) 
(region IV-b) C 2 > & C 3 > : P max = u\ . 

Finally, we compute u\ — u 2 , which is 

- vl = w + + - m ^ + Mh2 + it[2t + g) ) (5 ^ 15) 

where 

a 3 = h 6 + t(8t - g)h A - At 2 (g 2 + 5gt - 2t 2 )h 2 + 16#Y (5.16) 
(3 3 = h [h 4 + t(At - 3g)h 2 + Agt 2 (g - 2t)] . 
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One can show directly that u+ — v\ = when C 2 = 0. Also, it is straightforward to show 
that in g < 2t region v\ is always smaller than v\. Therefore, we can exclude v 2 + in regions 
VIII and IX. Combining all of these facts, one can express P ma x for 7 = n/2 case as follows: 

(i) g>2t (5.17) 
v\ C 2 > & C 3 > 

up remaining region 

(i) g<2t 
v\ C 2 > 

v\ C 2 < 0. 



p 

- 1 max 



P 

1 max 



Unlike 7 = case the whole parameter space is divided into the three applicable domains. 
Introducing the parameters u and v as Eq. fl4.13l) we plot the three applicable domains in 
the u-v plane in Fig. 4a. Around h = axis there are two domains, i.e. up and v\. Since vp 
and z/f go to fip and in the h — > limit, this guarantees that the h — > limit is consistent 
with same limit of 7 = case. The applicable domain for vp is little bit larger than the 
domain /j, P for 7 = case. The point (u = cos _1 (v^2/3), v = tan~ 1 (y / 3/2)) is shared by 
three domains. This point corresponds to 

\M = ||000> + ~ (|011) + |101) + |110» + z^|Hl>. (5.18) 

This is LU-equivalent with \W) = (1/^(1100) + |010) + |001>) as shown in Ref.Q- 

In Fig. 4b we plot the (u, v )-dependence of P ma x given in Eq. fl5.17p . Like Fig. lb the 
highly entangled states are represented as a valley in this figure. Fig. 4b seems to show 
that there exists an alley in the valley, which ends at u = v = tt/2. Along this alley so 
many highly entangled states are located. Comparing Fig. 4b with Fig. lb, one can realize 
that there are many more highly-ent angles states for 7 = it/2 case than 7 = case. This is 
mainly due to the fact that there are two LU-equivalent W-states when 7 = tt/2. 

Fig. 4c shows the geometric entanglement measure and the applicable domains simulta- 
neously in the u-v plane. Fig. 4c shows that around two W-states there are so many highly 
entangled states, which we would like to call W-neighbors. Especially, the neighbors of \4>w) 
in Eq. fl5.18l) gather along C3 = line. Besides the W-neighbors there are many highly en- 
tangled states around boundary of the applicable domains. These are the neighbors of the 



shared states 



22|, and we would like to call them the GHZ-neighbors. The GHZ-neighbors 
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0.0 0.5 1.0 1.5 

V 

FIG. 4: (Color online) Fig. 4(a) is a plot of the applicable domains for 7 = ir/2 case in (u, u)-plane. 
Unlike 7 = case there are three applicable domains in this case. The principal domain P max = v p 
is larger than P max = /-ip in 7 = case. This fact seems to indicate that the principal domain 
increases its territory with increasing 7. It is important to note that the domain P max = 1A is 
not reached to h = axis. This implies the consistency of the h — > limit. Fig. 4(b) is (u, in- 
dependence of P m ax- The highly entangled states forms a valley between two mountains. Fig. 4(c) 
is a plot of P m ax and the applicable domains in the (u, v )-plane. The boundaries of the domains are 
represented by black think line. Many highly-ent angles states reside around the boundaries and in 
the domain P ma x = v \. It is mainly due to the fact that there are two LU-equivalent W-states for 
7 = 7r/2 case. 

are slightly less-entangled compared to the W-neighbors. However, the number of the GHZ- 
neighbors are many more than that of the W-neighbors. 
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Finally, we consider the several special cases. First example is h = case. In this case 
C2 = —4:Q 2 t < and C3 = — (g — t) 2 (g + 2t) < 0, which results in identical expression 
with Eq. fl4.15p . Therefore, both results for 7 = and 7 = tt/2 cases coincide with each 
other in the h — > limit. Second example is t = case. It is easy to show that in 
this case P ma x — 9 2 when g > h and P m ax — h 2 when g < h. This is consistent with 
P nax (GHZ) = max(|a| 2 , \(3\ 2 ) when \GHZ) = «|000) + 

VI. EIGENVALUES AND GEOMETRIC MEASURE FOR 7 = vr/4 : NUMERICAL 
APPROACH 

In this section we will compute the eigenvalues and the geometric measure for 7 = 7r/4 
case. 

A. Eigenvalues 

For 7 = 7I-/4 Eq. OTOD reduces to 

2t(g + t) sin 6 cos tp + V2ht(l — cos 6) = A sin 6 cos if (6.1a) 

— 2t(g — t) sin 6 sin ip + \[2ht{\ — cos 9) — A sin 6 sin (p (6.1b) 
(g 2 - t 2 )(l + cos6) - h 2 (l - cos6) - V2htsm6(smip + cosy?) = Acos#. (6.1c) 
When 9 = 0, Eq. fl6.lap and Eq. (l6.1b|) are automatically solved and Eq.f l6.lcp gives 



\ = 2{g 2 -t 2 ). (6.2) 
Since s = (0, 0, 1) for this case, from Eq. (12.41) the corresponding eigenvalue is 

p 2 P = g 2 . (6.3) 



When sin^ 7^ 0, Eq. fl6.lap and Eq. (16. lb|) reduce to 



\-2gt-2t 2 X + 2gt-2t 2 , A . 

z = 7^7 cosv = 7^7 — sm ^ (6 ' 4) 

y 2ht v 2nt 

where z = tan(#/2). From Eq. (16.40 one can compute <p if A is known by using 

(A - 2t 2 ) - 2gt 

tanoj=— — . (6.5) 

Y {X-2t 2 ) + 2gt v ; 
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Deriving sin ip + cosy? from Eq. fl6.4p and inserting it into Eq. (l6.1c[) . one can derive the 
expression of z 2 in a form 

[(A ~ 2t 2 ) 2 - 4^ 2 ] (A - 2g* + 2t 2 ) 
" (A - 2h 2 )(\ - 2t 2 ) 2 - 8/i 2 t 2 (A - 2t 2 ) - 4^ 2 (A - 2h 2 ) ' K ' ' 

On the other hand, one can derive a different expression of z 2 directly from Eq. (I6.4j) 



z 



2 



(\-2gt-2t 2 ) 2 , 2 . , \(\-2t 2 ) 2 -4g 2 t 2 Y . N 

Equating Eq. (I6.6j) with Eq. (I6.7j) yields an equation for solely A: 

A/(A) = (6.8) 

where 

/(A) = A 4 - 2(h 2 + 4t 2 )A 3 - At 2 (2g 2 - h 2 - Qt 2 )X 2 (6.9) 
+8 [t\h 2 - 4t 2 ) + g 2 {3h 2 t 2 + 4t 4 )] A 
+ 16t 4 (g* - hg 2 h 2 - 2gH 2 - h 2 t 2 + t 4 ) . 

Eq. (16.81) guarantees the existence of the eigenvalue for A = as 7 = and 7 = tc/2 cases. 
In fact, one can show that there exists an eigenvalue corresponding to A = for arbitrary 
7. We have shown this fact in appendix A. 
When A = 0, Eq.([6J| and Eq.fJHZD reduce to 

2 (g 2 - 1 2 ) 2 g + t 

z = ttt~o — l?T tamp = -. (6.10) 

h 2 (g 2 + t 2 ) g-t 

Combining Eq. (l6.4p and Eq. (16. 101) . the possible solutions for 9 and ip are 

, 9 2 ~t 2 g-t . g + t 

z = ± coso9 = + — = sin oj = ± — =. (6.11) 

U^ff+t 2 ^ y/2(g* + 1 2 ) ^ ^2{g 2 + t 2 ) 

It is easy to show that both solutions in Eq. (l6.11|) gives a same eigenvalue, which is 

2 9 2 (9 2 + t 2 )h 2 + t 2 (g 2 -t 2 ) 2 
Po h 2 (g 2 + t 2 ) + (g 2 -t 2 ) 2 ' 1 ' ' 

Finally, let us consider /(A) = 0. It is worthwhile noting that at h — > limit /(A) = 

reduces to (A — 2gt — 2t 2 ) 2 (X + 2gt — 2t 2 ) 2 = 0. Therefore, the eigenvalues corresponding to 

/(A) = should coincide with /z 2 and /i 2 , for 7 = case, and with v\ and v\ for 7 = n/2 

case at the h — > limit. Equation /(A) = gives four solutions of A, say Ai, A 2 , A 3 and A 4 . 

We ordered the solutions by a fact that the h — *■ limit of Ai and A2 is —2t(g — t) and same 

limit of A 3 and A 4 is 2t(g + t). Then, the corresponding eigenvalues, say p 2 , p 2 ,, p|, and p 2 , 

can be computed numerically. 
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B. geometric measure 

FIG. 5: (Color online) Fig. 5(a) is a plot of the applicable domains for 7 = 7r/4 case. In this 
case there are two applicable domains. The principal domain P m ax 

= pp is little bit larger than 

Pmax = A*|> for 7 = and little bit smaller than P ma x = v p for 7 = it/ 2. This fact indicates that 
the principal domain increases its territory with increasing 7. Fig. 5(b) is (u, v )-dependence of 
Pmax- As 7 = case the highly-entangled states form a valley between two mountains. Fig. 5(c) is 
a plot of P m ax and the applicable domains in the (u, -u)-plane. Many highly-entangled states reside 
around boundary of the domains and near W-state. 

Using eigenvalues pp, pi derived analytically and pf {i = 1, 2, 3, 4) computed numerically, 
one can compute P ma x for the 7 = ir/4 case. Since each eigenvalue has its own available 
region, we checked this region by imposing Re [A] = 0, — 1 < sin6 < 1, — 1 < cos0 < 1, 
— 1 < siny? < 1, and — 1 < cosy? < 1. Although there are six different eigenvalues, the 
numerical calculation shows that only pp and p\ contribute to the geometric measure. This 
indicates that the whole parameter space is divided into two applicable domains. These 
two domains are represented m. u — v plane in Fig. 5a. The domains pp is slightly larger 
than domain pp and slightly smaller than domain up. This fact seems to indicate that the 
domain containing g = 1 extends its territory with increasing 7. 

Fig. 5b is a {u, v )-dependence of P ma x for 7 = 7r /4. Similarly with 7 = and 7r/2 cases, 
many highly entangled states reside at the valley between two mountains. Another highly 
entangled states reside around u — v — ir/2, which corresponds to W-state. The alley 
appeared in Fig. 4b does not appear in this case. This seems to be due to the fact that 
there is only one W-state in 7 = 7r/4 case. 

Fig. 5c is a (u, v )-dependence of P max and domains. As expected the highly entangled 
states are located around boundary and W-state. 

VII. CONCLUSION 

In this paper we have explored the effect of the phase factor in the geometric entanglement 
measure. We have chosen the most general three-qubit states which have symmetry under 
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the qubit-exchange. Our choice of the quantum states enables us to derive all eigenvalues 
and geometric measure analytically when the phase factor 7 is or n/2. It turns out that 
the 7 = 7r/2 case has three applicable domains while the 7 = case has two domains. Most 
highly entangled states reside around the boundaries of the domains and near W-state. 
Apart from the boundaries more and more the quantum states lose their entanglement and 
eventually, become the product states. 

Our result naturally gives rise to a question: what is a critical 7, say j c , which distinguish 
the two and three domains? In order to explore this question we have analyzed the 7 = 7r/4 
case numerically. Our numerical calculation shows that there are six different eigenvalues 
for 7 = 7r/4 case, but only two of them contribute to the geometric entanglement measure. 
Thus, there are two domains for 7 = 7r/4. 

We conjecture that emergence of the three applicable domains at 7 = n/2 is due to the 
two LU-equivalent W-states. In order to confirm our conjecture we checked numerically 
7 = 7r/3 and 7 = li7r/24 cases, which also give two applicable domains. We also checked 
the applicable domains for the partially symmetric quantum state 

|^) = ^lOOO) + f|011) + t|101) + f 3 |110) + e i7 /i|lll) (7.1) 

numerically when 7 = 0. This case also gives two applicable domains. Therefore, we 
conclude that the emergence of the three applicable domains is due to the appearance of 
additional W-state. 

In appendix we have shown that there exist eigenvalues for all 7, whose Lagrangian 
multiplier constant is zero. Although we conjecture that this is due to some symmetry of 
the quantum state we do not know the exact physical reason for the emergence of these 
solutions. It seems to be of interest to reveal the physical meaning of these solutions clearly. 
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Appendix A 

In this appendix we would like to show the existence of the eigenvalue /j/q, which corre- 
sponds to A = 0, at arbitrary 7. When A = 0, Eq. fl2.10j ) reduces to 

2/itcos7(l - cos 9) + 2t(g + t)sm6cosip = (A. la) 

2ht sin 7(1 - cos 9) -2t(g- t) sin 9 sin >p = (A. lb) 

(g 2 - t 2 )(l + cos 9) - h 2 (l - cos9) - 2htsm9cos((p - 7) = 0. (A.lc) 
The existence of /1q can be shown as following. First we derive 9 and ip by making use of 



Eq.l lA.lal) and Eq.l lA.lbl) . Then we show that the solutions 9 and <fi also solve Eq.l lA.lcll . 



Now, we consider only sin# 7^ case. Then from Eq.l lA.lal) and Eq. llA.lbj) it is easy to 
derive 

(g + t) sin 7 cos ip + (g — t) cos 7 sin (p = 0, (A. 2) 

which gives 

tany = — ^ tan 7. (A. 3) 

g-t 

Combining Eq. (1A.2l) and Eq. (1A.3l) . one can derive the solution for ip, which is 

jl 9-t . (g + t)tan-f 

coso9 = ± — , smo9 = =F — - (A. 4) 

V(g-t) 2 + (g + t)H a n 2 7 Y V(^-*) 2 + (^ + *) 2 t^ 2 7 



Inserting Eq. flA.41) into Eq. (lA.lbj) . one can derive sin# in a form 



sin 9 - t 2h (g 2 -t 2 )Vg 2 + ti-2gtcos2 1 

+ h 2 {g 2 + t 2 -2#tcos2 7 ) + {g 2 -t 2 ) 2 ' 1 ' } 

Inserting Eq. (|A.4f) and Eq. (jA.5p into the lhs of Eq. flA.lcl) . one can show straightforwardly 
that Eq. flA.lcj) is solved already by Eq. (1A.4I) and Eq. (IA.5j) . This guarantees the existence 
of 



In order to derive [1% explicitly we choose the upper sign in Eq. (1A.4I) and Eq. (1A.5I) . Then 
the components of the vector s becomes 

a singcosc -2% - t) 2 (^ + t) cos 7 

s x - sin cos - h2[{g2 + f2) _ 2gtcQsH + {g2 _ f2)2 (A.6a) 

• a ■ 2h(g-t)(g + t) 2 sm 7 

S * = Sm eSmif= h 2 [(g 2 + t 2 )-2gt C os2 1 ] + (g 2 -t 2 ) 2 (A - 6b) 

h 2 [(g 2 + t 2 )-2 9 tcos2 1 )-(g 2 -t 2 ) 2 
z h 2 [( g 2 + t 2 )-2gtcos2 1 ] + ( g 2 -t 2 ) 2 - 1 j 
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Inserting Eq. (lA.6|) into Eq. (I2.4|) and performing tedious calculation, one can show that /jLq, 

eigenvalue corresponding to A = 0, becomes 

2 = 9 2 h\g 2 + t 2 - 2gtcos27) + t 2 {g 2 - t 2 ) 2 
Mo h 2 (g 2 + t 2 -2gtcos2j) + (g 2 -t 2 ) 2 ' 1 ' 

It is straightforward to show that the choice of lower sign in Eq. (IA.4j) and Eq. (1A.5[) leads us 

to same expression of /x 2 ,. One can show easily that /i 2 , exactly coincides with /x 2 in Eq. fl3.lll) . 

v\ in Eq. fl3.29p and p\ in Eq. fl6.12p when 7 = 0, 7 = 7r/2 and 7 = 7r/4 respectively. 

Finally, making use of explicit expression of /x 2 ,, oen can derive the nearest product state 

\q)\q)\q') for //§, i.e. 



AB 



[q\(q\iI>) = t*oW) Ao{q\{q'\^) = fM>\q) Bc{q\{q'\ip) = Mo|«) (A.8) 



where j-0) is given in Eq. (12.31) . Since s is a Bloch vector of \q)(q\, one can show directly 
where 

2 22 g — t + ^ 

t =g +t — 2(?tcos27 cos 77 = — - — cos 7 sin 77 = — - — sin 7. (A. 10) 

Inserting Eq. ( 1A.9f) into Eq.( 1A.8l) it is straightforward to show that \q') becomes 

\q>) = 1 [{gh 2 £ 2 + t(g 2 - t 2 ) 2 e 2 ^} |0) + e^h(g 2 - t 2 ) {(g 2 - t 2 )e^ - 2tt) |1>] (A.ll) 

where A/" is a normalization constant, which makes \q') unit vector. 
For 7 = case the nearest product state becomes 

\v) = m = =~ (MO) -(g + Oil)) (A.12) 
+ (9 + *) 

|g') = 1 [{^ 2 + ffo + t) 2 } |0) + % 2 - t 2 )|l)] . 

^{gh 2 + t(g + t) 2 } 2 + h 2 (g 2 -t 2 ) 2 

It is interesting to note that (q\q') = when V>\ = 0, where T>\ is given in Eq. fl4.8p . 
For 7 = 7r/2 case \q) abd |g') becomes 

l?> = /u2 = „ 2 W) + % - (A.13) 

V h +(9- ty 

W) = 1 1 [{gh 2 - t( g - 1) 2 } |o> - z% 2 - t 2 )|i)] . 

{^/i 2 -t(^-t) 2 } 2 + /i 2 ((? 2 -t 2 ) 2 



It is interesting to note that (q\q f ) = when C3 = 0, where C3 is given in Eq. fl5.lip . 
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